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1. Introduction 

Rare decays of B or D mesons may play an important role in the discovery of new physics 
in the flavor sector. Transitions in which a b quark decays to an s quark proceed through a flavor 
changing neutral current. As these currents only occur at higher order within the Standard Model 
(SM), such decays are rare, and a small contribution from new physics beyond the Standard Model 
may be comparable to the SM contribution and hence observable. In this work, we focus on the 
B — )• Kl^l semileptonic decay process, which occurs via the transition of ft — )• sl^l^ at the quark 
level. References [|l|] and summarize both the experimental and theoretical status of B-meson 
decays. The BABAR Collaboration studied both B — )• Kl^l^ and B — )• K*l^l^ semileptonic decays 
in Ref. ||3|]. The Belle Collaboration published their results on the same process in Ref. [Q]. The 
CDF Collaboration studied the B K*}i+jX- decay in Ref. 

The B — Kl^l^ form factors are calculable from first principles using lattice QCD as there 
is only one hadron in the initial state and one in the final state. Recently, calculations using the 
MILC 2+1 flavor dynamical quark ensembles have been done for B — K*l^l^ by Liu et al. [||] and 
a study of B — /T* 7 form factors was done by Becirevic et al. [^. The Fermilab Lattice and MILC 
Collaborations presented some preliminary results for the B Kl^l^ form factors in Ref. [||]. 
Additional ensembles covering a wider range of lattice spacings have been analyzed since then. 
In this brief report, we show a more comprehensive (but still preliminary) analysis of the meson 
masses, form factors and chiral and continuum extrapolations. 



2. Theoretical Background 

An operator production expansion (OPE) analysis of B — )• Kl^l'^ shows that two currents, a 
vector current ly^b and a tensor current so^^qyb contribute to this process at lowest order [|]]. We 
study the vector current here and defer study of the tensor current to later work. The matrix element 
of the vector current can be expressed in terms of two form factors /+ and /o as: 

{K\rsfb\B) =f^ig2)(p^^+p^^-!±^gA+f,{g2y<^^,^ (2.1) 



where = Pg— p^- We study the form factors in the the B-meson rest frame, so only the kaon 
has non-zero momentum. The form factors f\\ and /l are defined as: 

/,=2^, (22) 

They are more convenient for our lattice calculation, and are related to /+ and /o by: 

/+ = [^1 + (ms - EK)fi] , (2.4) 

/o = [("^B - Ek) f II + (4 - rnDfj] . (2.5) 

mjj — mj. 

The form factor f±, as compared to gives the dominant contribution to /+, and hence to the 
experimental decay rate. 
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a (fm) 


amf^/amf^ 




Kb 


i. V TTl c n Tf* 


0.12 


0.02/0.05 


0.02/(0.0415, 0.05) 


0.0918 


2052 


0.12 


0.01/0.05 


0.01/(0.0415, 0.05) 


0.0901 


2259 


0.12 


0.007/0.05 


0.005/(0.0415, 0.05) 


0.0901 


2110 


0.12 


0.005/0.05 


0.005/(0.0415, 0.05) 


0.0901 


2099 


0.09 


0.0124/0.031 


0.0124/(0.0261, 0.0310) 


0.0982 


1996 


0.09 


0.0062/0.031 


0.0062/(0.0261, 0.0310) 


0.0979 


1931 


0.09 


0.00465/0.031 


0.0047/(0.0261, 0.0310) 


0.0997 


984 


0.09 


0.0031/0.031 


0.0031/(0.0261, 0.0310) 


0.0976 


1015 


0.09 


0.00155/0.031 


0.00155/(0.0261, 0.0310) 


0.0976 


791 


0.06 


0.0072/0.018 


(0.0072)/(0.0188) 


0.1048 


593 


0.06 


0.0018/0.018 


(0.0018)/(0.0188) 


0.1052 


827 



Table 1: Ensembles of QCD gauge field configurations used in the current B — )■ Kl^l^ work, amf^ and 
am^^^ denote the light and strange sea quark masses, am™' and am™' denote spectator and daughter quark 
masses in the b s transition. Both unitary and partially quenched kaon mass points are included in the 
chiral-continuum extrapolations. 



The lattice form factors and /j^ are computed numerically at several values of the lattice 
spacing and of the average up-down and strange quark masses. These results must then be ex- 
trapolated to the physical quark masses and the continuum using chiral perturbation theory (;^PT). 
SU (3) staggered chiral perturbation theory for heavy-light semileptonic form factors was suc- 
cessfully applied to the case of B — nlv decay in Ref. [[To|]. More recently, SU (2) ;^PT was applied 
to the extrapolation of D — nlv form factors on Nf = 2 lattices [^]. Some studies purport that 
SU{2) ;^PT may be a better effective theory for heavy-light physics projects [12]. Here we test 
both SU (3) and SU (2) formalisms in our B — ;> Kl^l^ chiral-continuum extrapolations. 



3. Numerical Simulation 

Our lattice calculations are done on MILC's Nf =2+1 flavor gauge configurations with asqtad 
improved quarks ||T3|]. The clover action with the Fermilab interpretation is used for the heavy 
quark [14]. MILC's ensembles cover many lattice spacings, light quark masses and volumes, which 
allows good control on the form factors' chiral and continuum extrapolations. The b quark mass is 
tuned close to its physical value as in Ref. [15]. In this report, we include results from the coarse 
(a ^ 0.12 fm), fine (a 0.09 fm), and superfine (a ?a 0.06 fm) ensembles. 



4. Numerical Results 

The first step of our analysis is to determine the B meson masses and kaon masses and energies 
from fits to the two-point correlators on every ensemble. States of both positive and negative parity 
contribute to the two-point correlators [|T0|]. We vary the number of states in the fits and select the 
fit range from fmin to ?max carefully to control the exited state contribution and obtain a good jc-value 
(confidence level) of the fit. Finally, the error on the mass is estimated via the standard jackknife 
method. 
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The second step of the analysis is to extract the form factors f\\ and f± from fits of the ratio of 
three-point and two-point functions. The three-point function is defined as: 

cf.;^(f,r;p^) = £e'''->x^^(o,o)y^(?j)^l(r,f)), (4.i) 

x,y 

where V^=isY^b and T is the location of the sink operator. Because we study the form factors in 
B-meson rest frame, only the kaon has non-zero momentum (pk)- In a finite volume, the kaon's 
momentum is discrete. We choose p = (0,0,0)27r/L, (l,0,0)27r/L, and (1, l,0)27r/L, where L 
is the box size. Higher momentum data are omitted due to their large statistical fluctuations. In 
addition, an iterative averaging trick is used to suppress the contribution of the states that have an 



alternating sign [ 10]. The ratio of two and three-point functions is defined as: 



'C2{t)C2{T -t) \ e t e ""B (T t) 

where C2 and C3 are averaged correlation functions [[Io|]. After multiplying R by the required 
renormalization constant, we obtain the continuum form factors:^ 



/|Pt^p^Zfz(/^f;^(^r), (4.3) 

/^«"t = p^^^^f;^(,,r). (4.4) 
Pit 



R3o(fpara) fit on a=0.12fm ml/ms=0. 01/0.05 lat. 



0.35 



0.3 



0.25 



0.2 



0.15 



0.1 



p-val=0.16, mom=000 
p-val=0.69, mom=100 
p-val=0.54, mom=110 



2 4 6 8 10 12 14 16 
t 



Figure 1: Example /|| fit on the coarse (a — 0.12 fm), flm//flms=0.01/0.05 ensemble. The y-axis is the ratio 
R-i O [c.f. Eq. ( j4-.l[ )] without any renormalization factors. 



'The calculation of the renormalization constants on the superfine ensembles has not been finished yet. We use the 
values from the coarse and fine ensembles to estimate these factors on the superfine ensembles. We also set p as 1 for 
this current preliminary analysis. 
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Figure |l| shows an example /ji fit on the coarse, ami /ams=0.0l/0.05 ensemble with p = 
(0,0,0)271 /L, (1, 0,0)271 /L, and (1, 1,0)271 /L. The y-axis is /?3o without multiplication by the 
renormalization constants. We fit R^fl with a constant term plus an exponential decay term. We 



choose the preferred fit range by fixing the size the fit interval, i.e., t„ 



n, but shifting the loca- 



tion of the minimum time slice ?min to obtain a stable central value and errors and a good /j-value. 
We plot the the result of the constant term and its error, which is consistent with the single plateau 



fit method used in Ref. | IC]. 



fpara NNLO SU(3) fit p=1e-23 



fperp NNLO SU(3) fit p=0.002 




0.02/0.05 
0.01/0.05 
0.007/0.05 
0.005/0.05 
0.0124/0.031 
0.0062/0.031 
0.0047/0.031 
0.0031/0.031 
0.00155/0.031 
0.0072/0.0188 
0.0018/0.0188 
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0.0018/0.0188 



Figure 2: /j| (left panel) and /i (right panel) chiral-continuum extrapolations with NNLO SU{'i) S;ifPT. 
Partially-quenched points are included in the fits, but are not shown in the figures for clarity. Open circles 
denote coarse data points, open squares denote fine data, and filled squares denote superfine data. Fit lines 
should pass through the data points of the corresponding color. 



fpara NLO SU(2) fit p=0.51 




0.02/0.05 
0.01/0.05 
0.007/0.05 
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0.0124/0.031 
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0.0018/0.0188 
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fperp NLO SU(2) fit p=0.51 




0.02/0.05 
0.01/0.05 
0.007/0.05 
0.005/0.05 
0.0124/0.031 
0.0062/0.031 
0.0047/0.031 
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0.0018/0.0188 
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Figure 3: /j| (left panel) and /x (right panel) chiral-continuum extrapolations with NLO SU{2) SjPT. 
Partially-quenched points are included in the fits, but are not shown in the figures for clarity. Open cir- 
cles denote coarse data points, open squares denote fine data, and filled squares denote superfine data. Fit 
lines should pass through the data points of the corresponding color. The cyan band shows the continuum- 
extrapolated form factor at the physical light-quark mass with statistical errors. 

To determine the form factors defined in the continuum with physical quark masses, a com- 
bined chiral and continuum extrapolations is applied. We use Staggered Chiral Perturbation Theory 
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(S;t^PT ) as the low energy effective theory. S;^PT accounts for taste symmetry breaking effects in 
the staggered quark action. NLO SU (3) S;t;PT supplemented by NNLO analytic terms was used 
successfully infi — )• tt andD — )• 7i decays [10, 16], but we find that it fails in the /ji fit for B — )• A'Z+Z^ 
process. Figure | shows the of resuh of the NNLO SU (3) S;tPT fit for /j| . The SU (3) fit gives a low 
p-value and incorrect behavior in the small Ek region, which may indicate the fact that the kaon 
is too heavy for SU (3) S;^PT . On the other hand, SU (3) ;^PT gives a reasonable extrapolation for 
f±, since the shape of f± is dominated by the B* pole and SU{3) S;^PT accounts for it correctly. 
However, the p-vahie. is still poor. 

Given these difficulties with SU (3) S;^PT, we switch to SU (2) S;^PT for the chiral-continuum 
extrapolations. Our SU{2) formula is inspired by the work of Ref. [pi]]. We take the heavy 
limit of the St/ (3) S;tPT result [§] to obtain the SU(2) reduction. Figure g shows that SU{2) SxPT 
works well with our data even at NLO and provides better control in the low Ek region. The results 
here are still preliminary. The detailed comparisons between SU (3) and SU (2) fits are still under 
investigation. 

After we obtain the continuum extrapolated /ji and f±, we can construct /+, which is crucial to 
the B — )• Kl^l^ differential decay rate. The plot of /+ is given in Fig. ^ Our simulation corresponds 
to the momentum transfer from 16GeV^ to 23GeV^. More study of the dependence of /+ in 
the full kinematic range using the z-expansion will be considered later [|T7|]. 



f+ vs. ot 



2.5 





1 - 



16 17 18 19 20 21 22 23 

(GeV^) 



Figure 4: Continuum extrapolated /+ for B Kl^l decay. The width of the band indicates the statistical 
error; in addition, we estimate that the systematic errors will be ^ 5%, depending somewhat on q^. 



5. Summary and Future Plans 

We report preliminary results from the study of the B — )• Kl^l^ semileptonic decay process 
with statistical error only. Current chiral-continuum extrapolations are done with NLO SU (2) 
S;tPT. More comparison of the SU (3) and SU (2) S;^PT will be studied. We will work on the 
systematic errors in the next step. These systematic errors come from the uncertainties in r\ (used 
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to set the scale) and physical quark masses, finite volume effect, chiral fits and so on. Finally, 
the S — ^ Kl^l~ process contains an additional form factor from the tensor current, which we will 
analyze in the future. 
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